First-Principles Study for the Anisotropy of Iron-based Superconductors toward 

Power and Device Applications 
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Performing the first-principles calculations, we investigate the anisotropy in the superconduct- 
ing state of iron-based superconductors to gain an insight into their potential applications. The 
anisotropy ratio 7a of the c-axis penetration depth to the ab-plane one is relatively small in BaFe2As2 
and LiFeAs, i.e., 7a ~ 3, indicating that the transport applications are promising in these supercon- 
ductors. On the other hand, in those having perovskite type blocking layers such as Sr2ScFeP03 
we find a very large value, 7a > 200, comparable to that in strongly anisotropic high-T c cuprate 
Bi2Sr2CaCu208-a. Thus, the intrinsic Josephson junction stacks are expected to be formed along 
the c-axis, and novel Josephson effects due to the multi-gap nature are also suggested in these 
superconductors. 

PACS numbers: 74.25. Jb,74.70.-b,71.15.Mb 



Since the discovery of iron-based superconductor 
LaFeAsOi-^F^ with T c = 26K by Kamihara et aZ.jJ, 
high-T c superconductivity has been reported in some of 
its family materials. All the superconductors in this fam- 
ily have FePn layers (Pn = P, As, or Se), in which the 
high-T c superconductivity emerges, and are classified into 
five groups according to the non-superconducting block- 
ing layers sandwiched between FePn layers as seen in 
Fig.HJ i.e., 1111 (e.g., LaFeAsO), 122 (e.g., BaFe 2 As 2 )[l, 
111 (e.g., LiFeAs) (jf, 11 (e.g., FeSe)[I], and the other type 
that contains very thick perovskite-based blocking layers 
(Sr 2 ScFeP0 3 )l. 

The quasi two-dimensionality of the electronic states 
due to the layer structure of the iron-based superconduc- 
tors is significant for their applications, since it causes 
anisotropy in superconducting properties such as the crit- 
ical current J c and the upper and lower critical fields, 
H C 2 and H c \. A convenient quantity characterizing the 
superconducting anisotropy is the anisotropy parameter 
7> defined as the ratio of the c-axis penetration depth 
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FeSe LiFeAs BaFe ? As 2 LaFeAsO Sr 2 ScFeP0 3 

FIG. 1: Crystal structures of typical iron-based superconduc- 
tors. 



to the in- plane one. From the values of 7a, one can 
extract important information about their potential ap- 
plications. For example, superconductors with small 7a 
are suitable for power applications using superconduct- 
ing wires, tapes or cables. On the other hand, those with 
extremely large 7a are expected to work as the intrinsic 
Josephson junctions [6jj|jjj which will be promising as a 
THz emission device 

It is well known that the penetration depth anisotropy 
7a directly leads to the anisotropy of H&, H c i and J c in 
single-band superconductors. In the multi-band systems 
such as the iron-based superconductors their relations 
become much more complicated and are still controver- 
sial floj |. But, in any case the parameter 7a describes 
the scale of the superconducting anisotropy. Thus, in 
this paper we perform a systematic first-principles study 
for 7a in the whole groups of iron-based superconduc- 
tors and elucidate the origin of the general trend in the 
superconducting anisotropy observed in various measure- 
ments such as torque, microwave surface impedance, J c , 

etc. (TIL E3, E3, E3, EE, US, E3, Though 

accurate values of 7a have not yet been settled experi- 
mentally, i.e., the reported values are still scattered, we 
find a clear trend in 7a among the five groups of the iron- 
based superconductors. Moreover, since we can predict 
the anisotropy ratio of a superconductor in which single 
crystals are not available, our method presented in this 
paper will provide a powerful tool for the exploration of 
the superconducting materials suitable for each applica- 
tion. Our calculations reveal that the anisotropy ratio 7a 
depends sensitively on the nature of the blocking layers 
and its values are widely distributed as seen in exper- 
iments, e.g., 7a ~ 3 in the superconductors belonging 
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TABLE I: Lattice constants a and c, the height hp n of As, P , 
or Se ions from the Fe-plane, and the anisotropy parameters of 
the resistivity and the penetration depth at zero temperature. 
a, c, and hp n are experimental values as the input parameters 
for the first-principles calculations, and j p and 7a are the 
calculated ones. 



i[A] c[A] fc Pn [A] 7p(0) 7a(0) 



FeSe 
LiFeAs 
BaFe 2 As 2 
LaFePO 
LaFeAsO 
Sr 2 ScFeP0 3 



3.7738 5.5248 1.4652 

3.7914 6.3639 1.6769 

3.9625 13.0168 1.3602 

3.9636 8.5122 1.1398 

4.020 8.7034 1.3238 



18.44 4.29 

9.06 3.01 

10.69 3.27 

17.34 4.16 

116.8 10.81 



4.016 15.543 1.1984 6.19 x 10 5 
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to the 111 and 122 groups and 7a ^ 200 in those with 
perovskite-based blocking layers. These results remind 
us of the variety in high-T c cuprates, e.g., 7a ~ 8 in 
YBa 2 Cu30 8 _5 and 7A > 100 in Bi 2 Sr 2 CaCu 2 8 _,5(Bi- 
2212) [22j. It is noted that the anisotropy in the 111 
and 122 groups is small in spite of the layer structure. 
The large value of 7a in the group with perovskite-based 
blocking layers strongly suggests the formation of the in- 
trinsic Josephson junction stacks along the c-axis. 

In this paper, we evaluate 7a for the typical iron- 
based superconductors from numerical results for the 
anisotropy ratio of the normal state resistivity at T = OK, 
7 P (0) obtained by means of the first-principles calcula- 
tions for the electronic states 23J, that is, the supercon- 
ducting anisotropy originating from the band-structure 
dependent Fermi velocities is considered. Let us first 
summarize the method employed in this paper to eval- 
uate the anisotropy parameters, 7 P and 7a- Assuming 
that the relaxation time of the conduction electrons is 
isotropic and independent of their velocities, we consider 
the quantity 7 P (0) = (v^fsA^fs to be equal to the 
ratio of the c-axis resistivity to the a-axis one. Here, 
v a and v c are the Fermi velocities, respectively, parallel 
and perpendicular to the FePn layers and (• ■ • )fs de- 
notes the average on the Fermi surface [2^]. The Fermi 
velocities are derived from the derivative of the band en- 
ergy, Vi = dek/dki, where the band energy e& is calcu- 
lated by the first-principles calculations. In our calcula- 
tions we use the first-principles density functional pack- 
age VASP 24j, which adopts GGA exchange-correlation 
energy 25 1 and PAW method (2(|. The electron self- 



consistent loops to obtain the charge density are repeated 
until the total energy difference becomes smaller than 
10 -6 eV. In these loops the spacing between nearest- 
neighbor fc-points is taken to be ~ 0.1-0.2 A . Once the 
charge density is determined, the energy bands are again 
calculated for more fc-points located with a much finer 
spacing, 79x79, in the xy-plane (the number of fc-points 
along the z-axis depends on the compounds) , in order to 
determine the Fermi surfaces and the Fermi velocities. 



(g) LaFeAsO 




FIG. 2: The Fermi surfaces of (b) Sr 2 ScFeP0 3 , (c) antifer- 
romagnetic BaFe 2 As2, (d) LaFePO, (e) FeSe, (f) BaFe 2 As 2 , 
(g) LaFeAsO, and (h) LiFeAs. The contrast on the surfaces 
indicates the magnitude of the Fermi velocity. The inset (a) 
shows typical symmetry points. 



The average (■ ■ ■ )fs is calculated by the standard tetra 
hedron method. In obtaining the anisotropy of the pen 
etration depth at zero temperature, we utilize the sim 
plest relation, 7a (0) 



V7p(0) [27|, which is convenient 
for systematic and comparative studies on the supercon- 
ducting anisotropy among various iron-based supercon- 
ductors, though its validity is limited in the multi-band 
systems 27J . The compounds that we examine in this 
paper are FeSe, LaFeAsO, LaFePO, BaFe2As2, LiFeAs, 
and Sr 2 ScFeP03. The calculations are done for tetrago- 
nal undoped compounds without any magnetic order, ex- 
cept for BaFe 2 As 2 . We calculate the electronic structures 
based on the observed crystalline structures, that is, the 
crystalline structures are not optimized, since it is known 
that the optimized structure deviates from the observed 
one in iron-based superconductors. 

The numerical results for the anisotropy parameters, 
7 P and 7a, in the iron-based superconductors employed in 
this paper are summarized in Table |TJ Let us discuss the 
1111-system, first. The experimental values in this group 
are widely distributed, depending on the compounds 

example, 7a ~ 15 — 20 



11, 13, 14 



distributed, depending 

EE EE El. I2E IliLFor . 



15j,l2l|, while 7a ~ 3 



in ReFeAsO (Re=Nd, Sm) [llj, 

PrFeAsO [2(|. Thus, we notice that there are two groups 
in the 1111-system, i.e., the one with relatively large 
anisotropy and the other with moderate anisotropy. In 
our calculations, we obtained 7a (0) = 10.81 for LaFeAsO, 
as shown in Table HI This value is close to that in 
the group with relatively large anisotropy. We also per- 
formed the calculations for stoichiometric PrFeAsO and 
obtained 7a (0) ~ 8.6, which is smaller than that in 
LaFeAsO as expected and also not far from the observed 
one in PrFeAsOi_ y . This difference between LaFeAsO 
and PrFeAsO can be mainly attributed to the difference 
of the main elements inside the blocking layers. On the 
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other hand, it is also noted that the calculated value in 
non-arsenic 1111-compound LaFcPO (7a (0) = 4.16) is 
much smaller than that of LaFeAsO, though the crystal 
structures are equivalent in both compounds. The Fermi 
surface shapes are different in these two compounds as 
seen in Figure [5] (d) and (g), that is, the Fermi sur- 
faces in LaFeAsO are almost cylindrical, while those in 
LaFePO are not, indicating that the two-dimensionality 
is weak in LaFcPO. From these results, one understands 
that the two-dimensionality differs among the 1111 com- 
pounds even if the blocking layers are equivalent. This 
remarkable feature comes from the fact that the band 
dispersion along the z-direction is sensitively related to 
the Pn's height from the Fe's square lattice plane. It is 
also noticed that the anisotropy of the 1111 compounds 
(7a si 10) is on the same order as in YB2CU3O8-5, that 
is, the transport application will be promising in the 
111 1-compounds. 

In the 122-compound BaFe2As2, we obtained 7a (0) = 
3.27. This value is smaller than that of LaFeAsO (1 Ill- 
system) and is consistent with the experimental one ~ 6 
[llj . which is generally smaller than that in the 1111- 
system. The weak anisotropy in the 122-compound is 
also understood from the feature of the Fermi surfaces. 
As seen in Fig 12 the Fermi surfaces in this compound are 
more winding along the z-direction compared with that 
in LaFeAsO, which leads to the weak anisotropy in the 
superconducting state of the 122-compounds. This indi- 
cate that the 122 system has high potential for transport 
applications. The 122-compound makes a tetragonal-to- 
orthorhombic transition, and a stripe-type antiferromag- 
netic (AFM) order appears in the low temperature phase. 
We also calculated the anisotropy parameter 7 P in the 
antiferromagnetic phase (fJ-Fe = 1-9/is) HH and found 
7 p (0) = 1.2, which is smaller than that of the tetrago- 
nal one being consistent with the experiments. In fact, 
the experimental value is 7 P ~ 2.5 [12j. Thus, one finds 
that the anisotropy becomes weaker in the AFM phase. 
Note that the reduction in the resistivity anisotropy orig- 
inates from the three-dimensional Fermi surfaces appear- 
ing in the AFM phase as seen in Fig. 2(c). These results 
clearly demonstrate the advantage of the present evalu- 
ation technique. In Table [I] we also list the anisotropy 
parameters of LiFcAs in the 111-systcm. The 7a value is 
comparable to BaFc2As2, that is, the Ill-system is also 
suitable for transport applications. This low anisotropy 
basically reflects the winding features of the Fermi sur- 
face like 122-systems as seen in Fig. 2(h). 

In the 1 1-compounds FcSe, which has the simplest 
crystal structure composed of a stack of only FePn lay- 
ers, the Fermi surfaces calculated in this compound shows 
clear curvature along the z-direction as seen in Fig. 2(e), 
indicating that the two-dimensionality is weak in this 
system. In fact, we obtained a relatively small value, 
7a(0) = 4.29, in the superconducting state. However, 
we note that this value is larger than those in the 122 
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FIG. 3: A schematic illustration of the intrinsic Josephson 
junction stacks in Sr2ScFePC>3. An alternate stack of FePn 
(superconducting) and blocking layers (insulating) are re- 
garded as a series of S-I-S Josephson junctions, ifi stands 
for the phase of the superconducting order-parameter of the 
ith electron energy band. The arrows represent possible tun- 
neling channels. 



and 111 systems and rather close to that in LaFePO. 
From this result one understands that metallic ions in the 
blocking layers enhance the three dimensionality in the 
122 and 111 systems. We emphasize that only the first- 
principles calculations can systematically derive such del- 
icate material-dependent difference in the anisotropy. 

Now, we discuss the anisotropy in the final group, i.e., 
the newly discovered iron-based superconductors with 
perovskite-type blocking layers, which are now under in- 
tensive exploration. The superconductors in this group 
do not show any magnetic order at all, while the su- 
perconducting transition temperature is rather high. In 
this paper we focus on a typical compound, Si^ScFePOs, 
which shows the highest T c at ambient pressure among 
non-arsenic iron-based superconductors. Since the block- 
ing layers in this compound are very thick as shown 
in Fig.l, strong anisotropy is expected in the normal 
state resistivity. In fact, a huge value 7 P ~ 10 6 is ob- 
tained, which is several hundred times larger than that 
in the other groups, i.e., 1111-, 122-, and 11- compounds. 
From this normal state value it follows, 7a (0) ~ 250, 
for the anisotropy in the London penetration depth of 
Sr2ScFeP03. Note that the value is comparable to that 
in Bi-2212 in high-T c cuprates. From this fact one may 
infer that the c-axis transport in this system is brought 
about by the electron tunneling between neighboring 
FePn layers, as in Bi-2212. In the superconducting state 
the Cooper-pair tunneling will be also expected to oc- 
cur between the superconducting FePn layers as in Bi- 
2212, that is, a single crystal of Si^ScFePOs may be re- 
garded stack of nano-scale Josephson junctions as 
schematically shown in FigUH which is called the intrin- 
sic Josephson junctions. In the intrinsic Josephson junc- 
tions the I—V characteristics under no external magnetic 
field show a remarkable feature, i.e., the multiple branch 
structure composed of many I — V curves @ , the num- 
ber of which corresponds to that of stacked intrinsic SIS 
junctions as shown in FigFJJ Hence, we naturally expect 
that the multiple-branch structure appears in the c-axis 
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FIG. 4: The penetration-depth anisotropy as a function of 
the distance between Fe-planes. 

I — V characteristics in the superconducting state in a 
single crystal of S^ScFcPOa. As for the damping nature 
in the Joscphson effects, this system should be under- 
damped one, because the paring symmetry is expected 
to be mainly full gapped s-wave one. It is also predicted 
that a new oscillation mode exists in the phase differences 
in addition to the Josephson plasma [34|. Moreover, the 
Josephson vortex has an internal structure (35j . These 
new features originates from the multi-tunneling chan- 
nels due to the multi-gap nature in these superconduc- 
tors as shown in Fig[3j The existence of the new phase 
oscillation modes will enrich the physics in the intrinsic 
Josephson junctions expected in the new compounds. We 
also mention that the superconducting Sr2ScFeP03 is a 
stoichiometric compound, that is, the superconductivity 
appears without doping. Then, it will be easier to make 
a stack of homogeneous intrinsic Josephson junctions in a 
large-scale, which is contrasted with the Bi-2212 intrinsic 
Josephson junctions, in which the doping for getting the 
high-T c superconductivity brings about disorders. 

Finally, we summarize our numerical results for 7a vs. 
distance between adjacent iron-planes in FigHl From 
this figure one finds a rough trend between the distance 
and the anisotropy, which indicates that the anisotropy 
increases with increasing the distance. However, our re- 
sults also suggests the existence of the other factors that 
affect the penetration depth anisotropy. As discussed in 
this paper, the anisotropy is sensitive to the distance be- 
tween the pnictogen ions and the Fe-plane, and also the 
elements forming the blocking layers. From these results 
we claim that the method presented in this paper, us- 
ing the first-principles calculations, is a powerful tool for 
evaluating the superconducting anisotropy. 

In summary, we investigated systematically the 
anisotropy of the London penetration depth for vari- 
ous iron-based superconductors on the basis of the first- 
principles calculations for their electronic states. It 
was shown that the anisotropy parameters, 7 P and 7a, 
are strongly dependent on the structure of the block- 
ing layers and pnictogen ions situated next to Fe lay- 
ers. Our numerical results well explain the variety of 
the 7a values observed in these superconductors. The 



122- and the Ill-compounds show weak anisotropy, being 
enough for transport power applications, while those with 
perovskite-type blocking layers have a very large value of 
7a , which is comparable to that in Bi-2212. We predict 
that the intrinsic Josephson junction stacks are realized 
in the superconducting state of these compounds. 

The authors wish to thank H. Ogino and J. Shimoyama 
for providing their recent experimental results and K. 
Tcrakura for illuminating discussion in first-principles 
calculations. The authors also thank N. Hayashi, Y. 
Nagai, M. Okumura, and N. Nakai for valuable discus- 
sion. The work was partially supported by Grant-in- 
Aid for Scientific Research on Priority Area "Physics 
of new quantum phases in superclean materials" (Grant 
No. 20029019) from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan. 



[1] Y. Kamihara et al, J. Am. Chem. Soc. 130, 3296 (2008). 
[2] M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 

101, 107006 (2008); K. Sasmal et al, Phys. Rev. Lett. 

101, 107007 (2008). 
[3] J.H. Tapp et al, Phys. Rev. B 78, 060505 (2008). 
[4] F.-C. Hsu et al, Proc. Natl. Acad. Sci. USA 105, 14262 

(2008) . 

[5] H. Ogino et al, Supercond. Sci. Technol. 22, 075008 

(2009) . 

[6] R. Kleiner et al, Phys. Rev. Lett. 68, 2394 (1992). 

[7] G. Oya et al, Jpn. J. Appl. Phys. 31, L829 (1992). 

[8] L. Ozyuzer et. al, Science 318, 1291 (2007). 

[9] K. Kadowaki et. al, Physica C 468, 634 (2008). 
[10] V.G. Kogan, Phys. Rev. Lett. 89, 237005 (2002). 
[11] R. Prozorov et al, Physica C 469, 582 (2009). 
[12] M. A. Tanatar et al, Phys. Rev. B 79, 134528 (2009). 
[13] D. Kubota et al. larXiv:0810.5623l 

[14] S. Weyeneth et al, J. Supercond. Nov. Magn. 22, 325 
(2009). 

[15] S. Weyeneth et al, J. Supercond. Nov. Magn. 22, 347 
(2009). 

[16] J. Jaroszynski et al, Phys. Rev. B 78, 174523 (2008). 
[17] M. Kano et al. |arXiv:0904.1418l 
[18] Z. Bukowski et al, Phys. Rev. B 79, 104521 (2009). 
[19] Y. Jia et al, Supercond. Sci. Technol. 21, 105018 (2008). 
[20] R. Okazaki et al, Phys. Rev. B 79,064520 (2009). 
[21] C. Martin et al, Phys. Rev. Lett. 102, 247002 (2009). 
[22] See e.g., J. C. Martinez at al, Phys. Rev. Lett. 69, 2276 

(1992), and references therein. 
[23] N. Hamada, T. Imai, H. Funashima, J. Phys.:Condens. 

Matter 19, 365221 (2007). 
[24] G. Kresse and J. Hafner, Phys. Rev. B 47, RC558 (1993); 

G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 

(1996). 

[25] J. P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 

77 (1996) 3865. 
[26] P. E. Blochl, Phys. Rev. B 50, 17953 (1994); G. Kresse 

and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
[27] V.G. Kogan, Phys. Rev. B66, 020509(R) (2002). 
[28] S. Margadonna et al, Chem. Commun. 2008, 5607 

(2008). 



5 



[29] J.H. Tapp et al, Phys. Rev. B 78,060505(R) (2008). 
[30] M. Rotter et al, Phys. Rev. B 78, 020503 (2008). 
[31] Y. Kamihara et al, J. Am. Chem. Soc. 128, 10012 
(2006). 

[32] C. de la Cruz et al, Nature 453, 899 (2008). 

[33] The calculated moment becomes larger than the observed 
one. This is a well-known problem in first-principles cal- 
culations for iron-based superconductors. We avoid fur- 



ther discussion on this discrepancy. See e.g., [3(| for its 
more details. 

[34] Y. Ota et al, Phys. Rev. Lett . 102, 237003 (2009). 
[35] Y. Ota et al, larXiv:0907.0277l 

[36] I. I. Mazin, M. D. Johannes, L. Boeri, K. Koepernik, and 
D. J. Singh, Phys. Rev. B 78, 085104 (2008). 



